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Abstract. The synchronization properties of two self-sustained quantum oscillators 
are studied in the Wigner representation. Instead of considering the quantum limit of 
the quantum van-der-Pol master equation we derive the quantum master equation 
directly from a suitable Hamiltonian. Moreover, the oscillators are coupled in 
incorporating an additional phase factor which shows up in the mutual correlations. 


1. Introduction 

In physics the ’synchronization’ stands for the emergence of temporally correlated 
behavior in composite systems. In the 20th century the main focus has been on the 
synchronization of periodically driven systems (see e.g. Ref. (^). In the past two 
decades, however, the search for synchronization has moved to chaotic systems |^. The 
interest in these systems stems from the fact that their temporal evolution sensitively 
depends on the initial condition, which is unfavorable for synchronization. Nevertheless 
it turns out that the trajectories of chaotic units can be synchronized under certain 
conditions. More recently the focus shifted again from classical to quantum systems. 
Here the presence of non-deterministic quantum fluctuations is expected to counteract 
synchronization, and therefore it would be interesting to explore under which conditions 
and in what sense quantum systems are able to synchronize. 

By dehnition synchronization is an irreversible phenomenon. Consequently, 
synchronization phenomena are only possible in systems with a broken time reversal 
symmetry. Since the standard quantum theory of closed systems is invariant under 
time reversal, it follows that quantum synchronization is only possible in open quantum 
systems that interact with the environment. On the other hand, interaction with the 
environment usually leads to decoherence of quantum effects j^, and this seems to be 
unfavorable for synchronization. 


The most common models for the study of decoherence in open quantum systems 
are the harmonic oscillator and two-level spin systems Setups where harmonic 
oscillators comprise both the system and the heat bath are studied under the name 
of quantum Brownian motion. It was found that the form of the coupling between 
system and bath determines the time scale of decoherence and the preferred states of 
the system (so called pointer states). On the other hand, spin-boson models describe 
situations where spin particles evolve under the influence of an harmonic oscillator bath. 
In this kind of models synchronization has been observed , depending on the coupling 
mechanisms between system and bath. The problem, however, is that the bath has a 
strongly dissipative influence: While two coupled oscillators will attune their frequencies 
over time, they will also both be damped due to the coupling to the bath. 

The observation that two non-driven quantum-mechanical systems cannot be 
synchronized is not too surprising. In fact, even in the classical realm synchronization 
is only possible if we are dealing with driven or self-propelled systems, simply because 
the desired synchronized state is an evolving one which depends on a continuous supply 
of energy to compensate dissipation. As we expect the same to be true in the quantum 
case, it is near at hand to look for synchronization in driven quantum systems out of 
equilibrium. 

The external drive could be implemented in different ways. For example, it could 
be introduced as a periodically changing driving force, but this would introduce an 
inherent time scale corresponding to the driving frequency which may interfere with 
other periodic properties of the system. For this reason we are interested here in freely 
evolving self-propelled quantum systems, which provide a mechanism to supply energy 
without imprinting a frequency. 

One of the simplest self-sustained oscillators in classical physics is the so-called van 
der Pol oscillator [^, which evolves according to the second-order differential equation 

Xe{l — x^)xX = 0 . (1) 

The drive is implemented in the second term, whose influence can be interpreted as a 
energy-dissipating friction for |T| > 1 and some kind of energy-injecting anti-friction for 
|x| < 1. As a result, the dynamics of the classical van der Pol oscillator converges into 
a limit cycle at |T| = 1 where both influences compensate one another. When suitably 
coupled, a ring of such classical oscillators is known to synchronize [^. It is therefore 
interesting to see what happens in the quantum case. 

Recently a quantum analog of the van der Pol oscillator has been studied in 
Refs. [^-[^. For example, in Ref. it was found that with increasing coupling 
between two oscillators the weight of the joint-probability distribution P{xi,X 2 ) becomes 
gradually concentrated near the line Xi = X 2 , indicating phase correlations. In 
conclusion the authors remark that strict frequency locking (as apparent in classical 
systems) cannot be found in quantum systems due to quantum noise. However they 
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state that they find strong frequency entrainment. 

In this paper we extend these previous works. First we will review the master 
equation for the quantum van der Pol oscillator, presenting a different derivation. 
Having studied the time evolution of a single unit we then focus on two coupled van 
der Pol oscillators in the second part. In particular, we introduce a generalized coupling 
involving an additional phase factor and calculate the steady state solution. As usual, we 
visualize the synchronization properties by using Wigner functions. Moreover, we show 
that the visualized results are in agreement with a direct evaluation of the correlation 
between the oscillators. Finally we calculate the entanglement of the two oscillators. 


2. Setup of the master equation 

Let us first consider a single quantum van der Pol oscillator. Like an ordinary harmonic 
oscillator this system possesses a discrete energy spectrum. Because of the coupling to 
a heat bath, the energy levels are randomly populated. This can be described in terms 
of a time-dependent density matrix p{t). For high energies, where the distribution 
over high energy levels is quasi-continuous, we expect the system to behave classically. 
Conversely, quantum effects are expected to play a role in the limit of low energies. 
Following Refs. (^-[^ we will consider here the quantum limit at very low energies. In 
this limit the physics is dominated by the ground state and the first excited state while 
all higher energy levels can be neglected. This effectively leaves us with a two level 
system. 

The time evolution of a van der Pol oscillator in this quantum limit is given by the 
master equation 

pit) = -i [H, p{t)] + K {2a+p{t)a- - p(t)}) (2) 

+ A (2cr“p(t)cr+ - {cr+cr“, p(t)}) , 

In this equation and a~ are the creation and annihilation operators acting on the 
two lowest levels of the oscillator. The free system Hamiltonian is given by 


H = uja~^a , (3) 

where cu is the oscillator frequency. Moreover, k and A are the coupling constants of the 
interaction with the bath. 

This equation can be derived in two different ways. Following one can start 
from the general quantum master equation for van der Pol oscillators, p = —i [H, p] + 
KiC [al] P + K 2 C [a^] p, where C denotes the Lindblad superoperator. By taking the limit 
^2 —>■ 00 one arrives at Eq. 0 with K = Ki and A = 2ki. 

Alternatively one can start with the Hamiltonian of a spin system in a spin bath 
similar to the approach in Ref. [^. The time evolution is then obtained by assuming 
that the system interacts with one bath spin at a time for a short interaction period. In 
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this approach the two terms in Eq. Q with k and A are generated by bath spins that are 
initially in the ground or excited state, respectively, before the interaction takes place. 
In this case we find that the rate of the incoming particles in the ground and excited 
state can be chosen individually. This means that the ratio k,/\ does not necessarily 
need to be 1/2. Moreover we see that k and A can be interpreted as the probability 
of finding a bath spin in the ground or excited state. The ratio of these parameters 
can therefor be translated into a temperature of the bath by a Maxwell-Boltzmann 
factor 
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A 

K 


exp 


huj 


( 4 ) 


In physical terms any bath temperature is allowed in this second approach as opposed 
to one specific temperature in the first approach. Technical details of the derivation can 
be found in [Appendix A[ 

To study correlations and synchronization properties of quantum systems one needs 
at least two self-sustained oscillators which are mutually coupled. The corresponding 
master equation can be obtained by adding the contributions of both oscillators 
supplemented by a suitable coupling term: 
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p{t) = -i[H,p{t)] +kY] {2a+p{t)a;- - , p{t)}) 


i=l 


( 5 ) 


+ A (2a^ p{t)a+ - {a+a^ ,p(t)}) 

i=l 

+ V (2cp{t)c^ — c^cp{t) — p(t)c^c) 


Here the scalar V determines the strength of the coupling between the two oscillators. 
In the following we choose A = 2k, E > 0 and u is the same for both. For the operator 
c occurring in the coupling term we use the choice c = 
in detail in section [5l 


(72 , as will be discussed 


3. Phase Space Representation using Wigner functions 


Quantum systems described by the master equation ([^ can be realized in different 
physical setups. One practical and experimentally advanced realization are quantum 
optical experiments, where the qubits are represented by single photons. In this research 
area the Wigner quasi probability function is a common tool to visualize the qubit’s 
quantum state p and its evolution. It also provides a convenient way to compare the 


quantum case with the classical behavior. Following 11 we now derive a suitable form 
of the Wigner function. 
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First we define the quadratures 

cr++cr" ^ .cr+-cr" 

Since a~ = -^{q + ip) it is clear that a~ encodes the complex amplitude of the oscillator 
in phase space and that the quadratures represent its real and imaginary part. For 
an optical realization the quadratures correspond to the in-phase and out-of-phase 
component of the electric field amplitude. Since they can be regarded as position and 
momentum of an oscillator, it is therefore reasonable to represent them in a classical 
phase space. 

For a pure quantum state \'ip) the Wigner function is given by 

/ OO 

ip*{x + y)'tp{x — y) dx. (7) 

-OO 

where 'ip{x) = {x\'p) and q\x) = x\x) and consequently 

ip*{x + y)'p{x -y) = {x- y\ip){pj\x + y) = {x- y\p\x + y). 


Realizing that the pure state p = \'p){'p\ is mapped linearly to the Wigner function W, 
it is clear that the same relation will hold for mixed states as well. Writing the mixed 
density matrix in the orthonormal energy eigenbasis p = '^rn n Pm,n\'^) Wigner 

function can be expressed as 

bF(x,p) Pm,nlRm,nP); (8) 

m^n 


with 


Wm,n{x,p) = I {x- y\m){n\x + y) e^*™dy 
Pm^x - p)0* (a: -F y) dy. 


^ —OO 

poo 


( 9 ) 


where (pnix) = {x\n). Using the common representation of the Heisenberg algebra 
(g —7- x and p these wave functions are just the harmonic oscillator 

eigenfunctions 


Pnix) 


1 


7r4”(n!)^ 


Hni^X^ , 


( 10 ) 


where Hn{x) are the Hermite polynomials. This allows us to compute the first four 
components of the Wigner function 

Wop{x,p) = 

Wq^i{x,p) = -\- ip) (11) 

hFyo(a^,p) = — ip) 

Hi,i(a:.p) = -|- 2p^ — 1). 
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Re [ Wi,o [x, p] ] 


Im [ Wi,o [x, p] ] 



Figure 1: 
Eq. (P. 


Components Wm,n{x,p) of the Wigner function of a two-level system in 


which are plotted in Figj^ As one can see, Wo,o and Wi^i are rotationally invariant in 
phase space while for the other two functions this symmetry is broken. 
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Figure 2: Wigner probability density of a single quantum van der Pol oscillator at 
different times. The initial condition is given by equation (13). For long times t the 
distribution approaches the steady state solution. 


4. Properties of a single quantum van der Pol oscillator 


Let us first consider a single quantum van der Pol oscillator. Since the master 
equation i§ can be interpreted as a set of three linear differential equations it can 
be solved exactly. For t ^ oo one finds that the density matrix p{t) always approaches 
the steady state solution 

Pstat = = 2|0)(0| + -|1)(1|. (12) 


This state is a probabilistic mixture of the ground state and the first excited state. 
Because of (12) the corresponding Wigner representation is symmetric under rotation 
in phase space. 


Turning to the time evolution it would be interesting to see the actual oscillation 
of the quantum oscillator in the Wigner representation. To this end we have to start 
with a probability distribution which is peaked at a certain point in phase space away 
from the origin. For example, starting with the initial (pure) state 


P(0) = 5(|0>-i|l>)((0|+*(l|) 


(13) 


which is not symmetric under rotation in phase space, we obtain the time evolution 
shown in Fig. As can be seen, the maximum of the probability density rotates 
counter clockwise and moves to the center until the steady state distribution (12) is 
reached. 


Why does the dynamics evolve into a rotation-symmetric probabilistic combination 
of |0)(0| and |1)(1|? This can be explained as follows. As one can see from equation 0 
the bath is coupled to the system via the creation and annihilation operators and 
a ~. These operators switch between the two levels and destroy any quantum-mechanical 
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superposition of |0) and |1). Therefore the phase information is destroyed by decoherence 
and the system is localized in a probabilistic combination of the ground and excited state. 
The mixing ratio (1:2) between the two states is given by the ratio of k and A in Eq. ([^. 


5. Coupled oscillators 


Let us now study possible synchronization effects between two coupled quantum 
oscillators, as described by Eq. ([^. In the literature the operator c is usually defined 
as c = ai — a 2 [^. In the present work, however, we use the more general form 

c = erf + e*'^crf (14) 


with an additional conpling phase parameter 0. The conpling operator c can be 
interpreted as a simultaneous transition between ground and excited state of both 
oscillators which is indneed by the heat bath. Depending on the coupling phase 0 
one of the oscillators acquires an additional phase factor relative to the other. As we 
will see below, this phase will show up in the phase probability distribution of the two 
oscillators in the steady state. 

Again it is straight-forward to compute the steady state. It depends only on the 
ratio 7 ^ while it turns out to be independent of oj: 




P 


7+3 


7(87+27)+27 

0 


0 

0 


0 

(7+2) (7+3) 

7(87+27)+27 

e~^‘^7(7+l) 

7(87+27)+27 

0 


0 

__d^7(7±l)_ 

7(87+27)+27 

(7+2)(7+3) 

7(87+27)+27 

0 


0 

0 

0 


\ 


2(7(37+8)+6) 
7(87+27)+27 / 


(15) 



(a) Probability density of the phase 
difference AO for given coupling 
phase (j). 


<p1 0P2 > 



(b) Expectation value of the momentum correlation function 
given in Eq. 0 for K = 1. Initially there is an 

approximately linear increase. 


Figure 3: Different measures of synchronization for two coupled oscillators. 
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5.1. Influence of the coupling phase 


The Wigner function can be defined in the same way for two oscillators. However, 
instead of using Xi,X 2 ,Pi, and p 2 as arguments we make a transformation to polar 
coordinates with 


= V2(xj + ipj) (i = 1 , 2 ) 


(16) 


as described in the Appendix of Ref. j^. Using this polar representation we can ignore 
the magnitude and only look at the probability distribution for the phase difference 
A6 = $2 — 0i of the two oscillators. 

We find that the probability is maximal for A9 = (p — ir. This can be seen in Fig. 
where the probability density of the synchronization phase A9 is plotted in dependence 
of the coupling phase f for fixed values U = 10 and k = 0 . 1 . 

Synchronization can also be detected by measuring the momentum correlation 
fnnction which is defined as the expectation value 


(Pi ® P 2 ) = Tr[(pi ® p 2 )p] 


7(7 + 1 ) cos( 0 ) 
7(87 + 27) + 27 ■ 


(17) 


This correlator depends on the coupling phase f and vanishes if the oscillators are 
coupled out of phase (meaning that 0 = |). Moreover, as can be seen in Fig. 3b and 
Eq. (17) it is bound by \ {pi ^ P 2 ) \ A Obviously, the behavior of this correlator is in 
agreement with the previous findings obtained by the Wigner function. 


5.2. Entanglement 

Ef[V] 



Figure 4: Entanglement of formation for two coupled oscillators (steady state solution 
of eqnation @) plotted over conpling strength V. 
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An important property of coupled quantum mechanical systems is the emergence 
of entanglement. For two two-level systems (i.e. qubits) the common entanglement 
measure is the entanglement of formation . It is defined by the formulas 


Ef{p) = S 




with 5'(a:) = —a:log 2 a: — (1 — a:) log 2 (l — a:) (18) 

and C{p) = max(0, Ai — A 2 — A 3 — A 4 ), 


where {Ai, A 2 , A 3 , A 4 } are the eigenvalues of the operator p((T^ ® 
by magnitude in descending order. In the present case we have 

^ 4(7(7+ 3 ) + 3 )^ 

^ (7(87 + 27 ) + 27)2 

, _ 2(7+ 3 )(7(37+ 8)+ 6) 

“ (7(87 + 27 ) + 27)2 

4(27 + 3)2 

" ( 7 ( 87 + 27)+ 27)2 


)p*(cr^ ® (j^) sorted 


(19) 

( 20 ) 
( 21 ) 


SO that Ai — A 2 — A 3 — A 4 = • Inserting this expression into the formulas 

listed above, the entanglement of formation can be computed directly. The result is 
shown Figure As can be seen, the entanglement vanishes for small values of 7 . The 
onset of entanglement happens at a coupling ratio of 7 = 2^/3 ~ 2.847. For large 

values of 7 the entanglement reaches an upper limit of Ef{p) = Afl + ~ 0.01118. 

This means that even for very large coupling strength the entanglement is weak. A 
maximally entangled two-qubit state, e.g. a Bell state, has Ef = 1. 


5. 3. Eigenvalue spectrum of the reduced density matrices 


In order to study the influence of mutual coupling on a single oscillator, we compute 
the reduced density matrix: 


( 7V^+2iyK+18+ 
8V2+27yK+27K2 

0 


0 

(U+3k)^ 

8V2+27yK+27K2 


( 22 ) 


In the limit —>■ 0 this reproduces our single oscillator result (compare section . For 
large values of V the eigenvalues of the single oscillator tend to | and 4, as shown in 
Fig. [5} The parameter k defines how quickly this transition happens. 


6. Conclusions 

In this work we have studied the synchronization properties of two self-sustained 
quantum oscillators coupled by an interaction term involving a phase shift. We confirm 


10 













Figure 5: Spectrum of the partial density matrix (equation (22)). The horizontal lines 
are the eigenvalues of the uncoupled oscillator. 


that the interaction leads to phase correlations between the systems. In particular, these 
correlations depend on the phase shift in the interaction term. These correlations have 
been visualized by using Wigner functions and evaluated explicitly. The correlations are 
a property of the stationary state, meaning that they are not destroyed by decoherence. 
Thus the system shows a stable correlated behavior although it does not synchronize 
perfectly in the classical sense. We have also analyzed the entanglement of the oscillators 
as a function of the coupling strength in the stationary state. While correlations are 
present for any value of > 0 we find that entanglement only sets in at a well-define 
threshold of the coupling strength. 

It would be interesting to extend this analysis to more than two oscillators. For 
example, in the classical case it is known that four oscillators synchronize partially in 
subgroups, and it would be interesting to find out whether this feature survives in the 
quantum limit. Moreover, it would be interesting to understand the meaning of the 
eigenvalues of the reduced density matrix. 


Appendix A. Alternative derivation from a spin system in a spin bath 

Eq. ([^ can also be obtained in a different way by following Englert and Morigi , who 
derived a master equation very similar to the quantum van der Pol master equation, 
however, lacking the square in the second Lindblad term. 

If we perform an analogous derivation with an oscillator Hamiltonian of a spin 
system in a spin bath it gives us right away the master equation in the quantum limit ([^ . 
Starting point of the derivation a single spin representing the system s and an arbitrarily 
large number of spins representing the environment e. 
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For the interaction part we make the following assumption. We can imagine that 
at different times individual bath atoms touch the system and interact with it. Then 
the system evolves freely again until the next bath spin arrives. We assume that only 
one of the environment spins interacts with the system spin at any time. Therefore, we 
consider the Hamiltonian 


H = uj 




) + 9 {(■ 




H 


free 


Hi, 


(A.l) 


where only one of the environment spins appears in the interaction part. 

In the end, when we put together the master equation, we will account for the 
different bath spins that interact with the bath at different times by performing an 
appropriate average. 

In the literature on quantum optics a similar Hamiltonian is studied under the name 
“resonant Jaynes-Cummings interaction in the rotating-wave approximation”, where 
the system is modeled by a harmonic oscillator instead of a spin. The first term Hfree 
describes the free evolution of the two spins. For simplicity we assume that they both 
have the same frequency uj. The coupling between the two spins is described by the 
second (exchange) term Hint controlled by the coupling strength g. The time evolution 
of the density matrix p is given by von Neumann’s equation 

= i [ft, H]. (A.2) 


Let us now look at a timespan [t, t + r] during which the interaction between the system 
and a single bath electron takes place. For this we first note that the two terms in the 
Hamiltonian (A.l) Hfree and Hint commute: 

= [(U ((j!cr„ 


[Hf reo H, 


inti 


alcT^ 








)] 


(A.3) 


= ujg{[a\a^,ala^] + [ala^, ala^] + + [ala^,ala^]) 

= <^9 {<^1 [(^ s , ^l] + crl [^ 1 , ^ e ] 

= ^9 = 0 • 


Under this condition the state of the total system at time t + r reads 


pt+r = (A.4) 

— g freef 
= ree r pte^'^‘^~\ reeT 

where we introduced the notations 7 = + crjug) and f = gr. We now rewrite the 

term which amounts for the change introduced by the coupling. The first 

exponential can be written as 


int'T’ j 


^iHfr 


g #7 _ cos((/) 7 ) + i sin((/) 7 ) = cos 



sin 
+ i'y — 



\/V 


(A.5) 
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It is convenient to nse this form because we can apply some general rules, namely 



= 1 

(A. 6 ) 


= 0 

(A.7) 

(t/((tV) 

= 

(A. 8 ) 

crVpV) 

= t^y(o) 

(A.9) 


= 

(A.IO) 

cr/(crcr'^) 


(ATI) 


These rules apply to functions / of operators in general. Applying the basic relations 


of fermionic operators (A. 6 ) and (A.7) we can calculate 7 ^ = + a^a^a^a^. Let 


us use the rest of the rules to transform a function into a suitable form, which we 

can then employ for the next step of the derivation. 

^( 7 ^) = {ala^ + (jgfj]) F{ala^(T^(Tl + a^alala^) (A.12) 


= 1 


= ^(1 ® + 0 • + a^al F{0 ■ a^al + 1 ® ct]ctJ 

= F{a,al) + a^al 

7 ^( 7 ^) = F{ala^a^al + a^alala^) 

= 4 (^s + Fi(^s4) 


(A.13) 


We can apply this to (A.5): 


sin 


cos 


cr^CTg cos 




(A.14) 


■e4 


crlcF^ 


c 


c 


sin! (l)^J a ^al 


+ ia^ al- 


sm 


+i(jl 


a^al 


a.a's 


<y,a's 


(J. 


S' 


St 


In the following we use the abbreviations ((7, (7, and S) introduced in this 


equation (A.14). 


Up to here our derivation depends only on algebraic relations and is independent 
of the choice of the density matrix. To be able to proceed with the derivation, we now 
choose a representation for the (jg-operators. 

Consider a vector representation of the bath atom where |0) = (0,1)^ and 
|1) = (1,0)^. Using this representation the cTg-operators are actually 2 x 2-matrices. 
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With this notation (A. 14) simplifies to 

sin(^0V^) 


= cos 


T + ^ 7 - 



(A.15) 


where the final matrix is actually a 4 x 4 matrix. 

Next we make some assumptions about the density matrix. Before the interaction 
(at time t), the system and bath particle are independent. In mathematical terms this 
means that they are separable and the total density matrix can be written as a tensor 
product of the two contributing terms pt = 0 Pt ^^. Ultimately we are only interested 
in the reduced density matrix of the system. Tracing out the bath particle we obtain 


47 = Tr, [. 


-iHr 


Pte 


iHr 


(A.16) 


— lOJCr'sCF 




Tr, 




^lUJ(Ts(JgT 




c 

-iS 


-iS^ 

c 







At this point we recall our assumption about the bath spin that we made in the 
beginning, that it is one of many spin particles in the bath. Now we specify the state 
of this bath spin at the time when the interaction starts. It can either be in the g round 

state p^p ~ ^ 0 1 ^ excited state p^p ~ 0 0 ^ time evolution (A.16) 

yields a different result for each of these states. Therefore we do the next steps for both 
cases in parallel. 


For a particle arriving in the ground state we have 

0 0 

0 4 " 

and for the excited state we have 

C -iS^ 

-iS C I \ 0 





s'p\’’s -iS'p-pC 

iCpf'S CpPc 

Cp'-pC iCp'pS' 
-iSpi*>C SpfSt 


We can finally calculate the trace by inserting this into the time evolution ( |A.16[ ) for 
the ground state 

icp[^^s cp^:^c 


pit = e-*‘^"'"«"Tre 


^iuja'sCFgT 


(A.17) 


= e 


—lUOCFsCf^T 




^lUJa'sCFgT 


and the excited state 


— ^-i^crlagT 

Pt+T — ^ 


Tr. 


c4''c 

-iSp'pC 


iCpPS* 

spY's' 




(A,18) 
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_ ^-iujcr^cr^T 


Cp\"^C + Sp\"^S^ 




To proceed further we expand these terms around 0 up to 2nd order in (/). This makes 
our derivation only valid for weak interactions. We pick up the derivation in (A. 17) and 
explicitly write out the terms C, C, S and ^ before we expand the sine and cosine 
terms in 6. 


(s) _ 


Pt+T ~ ^ 


lUJa's CTgT 


= ^-i^cr'sCr^T 


^-lUd'sCFsT 


^-lUJcr'sCr^T 


S^p^t'^S + Cp^t^C 


^lUJCr'sCFsT 


(A.19) 


smi (pxja^a's 




(s) t 
(^sPt 


sini (pxja^a's 


a^al 


+ cos ((j)\J alcr^ j cos i (j)^Jala 


^iuja^a^T 




^iua'sU^T 


1 


2 . V 2a,pi ’al - ala^p] - p\ V](j,) + p\ 


is) 


^iiXia'sa^T ^ 


Applying the same process for the interaction with bath particles initially in the excited 
state we get a time evolution 


(*) _ 


Pi+ 


luja's OgT 


Cp^^C + 


^lUJCr'sCFsT 


(A.20) 


= e 


—iLuaLa^T 


-y (2alpi"^a^ - - pJ^V^crJ) + pi 


( 5 ) 


lUjala^T 


Note that the 0th order term in both equations is the free evolution of without the 
bath. The change induced by the bath alone is 

- ala^pi"'^ - (A.21) 

or 

(2alpi"'>a^ - (T^alp["'> - . (A.22) 

We can further assume that the arrival times of these bath spins are independent, this 
means they follow Poissonian statistics. The average time between the arrival of two 
bath spins (i.e. the inverse of the arrival rate) is assumed to be much larger than the 
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interaction time r. This is necessary to ensure that not two bath spins interact with 
the system at the same time. If we specify a n arr ival rate k. for ground state atoms and 
A for excited state atoms we can now write (A.2) for : 
d 


dt 


(s) 

Pt 


TT 

Pt 5 ^ 


(A.23) 


(s) f 

Pt 


+ \ex 


^^sPt - (^WsPt - Pt 

o t (*) t (*) (*) t 

^<^sPt - (^s<^sPt - Pt 


This is the same as Eq. with just a little difference in notation: Here and a are 
used instead of and a~. 
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